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Abstract 



< 

In the preceding articles we considered fractional integral transforms involving one real scalar variable, 
one real matrix variable and real scalar multivariable case. In the present paper we consider the multivariable 
case when the arbitrary function is a real- valued scalar function of many pxp real matrix variables X\, X^. 
Extension of all standard fractional integral operators to the cases of many matrix variables is considered, 
along with interesting special cases and generalized matrix transforms. 

1. Introduction 

In all the preceding papers in this sequence the basic claim is that fractional integral operators are of two 
qq ■ kinds, the first kind and the second kind. The first kind of operators belong to the class of Mellin convolution 

| of a ratio and the second kind of operators belong to the class of Mellin convolution of a product. We will 

give the following formal definition of fractional integral operators of the first kind and second kind. The 
| following standard notations will be used in this article. All the matrices appearing are pxp real symmetric, 

and further, positive definite (denoted by X > O) unless otherwise specified. \X\ denotes the determinant 
C*~) ' of the matrix X, tr(A) denotes the trace of X, 



dX = J^J /\dxij for a general X 



= J^J Adxij = J^J Adxij when X is symmetric, X = X' 

i>j j>i 

that is, the wedge product of all differentials in X. Jo<a<x<b f(X)<iX = Jq f(X)dX will mean the integral 
over all positive definite X, such that A>0,X>0,X-A>0,B>0,B-X>0 where A and B are 
constant matrices. In this notation O < X < I will indicate that all eigenvalues of X are in the open interval 
(0, 1). All functions of matrix argument considered in this paper are real- valued scalar functions, whether 
the argument is one matrix or more matrices. The real matrix-variate gamma function will be denoted and 
defined as 

V p (a) = n PjE ^T(a)T(a - ~)...r(a - »(a) > (1.1) 

Also, will denote the positive definite square root of a real positive definite pxp matrix X . A prime 
will denote the transpose, that is, X — X' means X is symmetric. The following standard Jacobians will be 
used frequently in this article. For more details and for more Jacobians see Mathai (1997). 
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Y = AX A', \A\ + 0, X = X' => dY = \A\ p+1 dX (1.2) 
_ ! _ J |r|" 2p dy for a general X , . 

dA ~ \ |y|-(H-D d r for x = x' (1 - d > 

Definition 1.1. Fractional integral operators of the first kind in one scalar or matrix variable case A 
fractional integral of the first kind of order a and of one scalar or matrix variable is a Mellin convolution of 
a ratio with the first function fi(xi) is of the form 

/i(»i) = ^ (ai)( i- iCl)a " 1 ,8i(a)>0 (1.4) 
r(a) 

in the real scalar variable case, 

■KxQiJ-A^r^ p-i 

in the single matrix variable case, where 4>\{xi) or 0(-X"i) is a specified function, and 

f 2(^x2) — (j>2 (#2) f{%2) where <j>2(x2) is a specified function and f{x2) is an arbitrary function, in the 
real scalar case and 72(^2) = faiX^) f (X^) in the real matrix case where X2 is a p x p real positive definite 
matrix, so that the fractional integral of the first kind of order a is given by the Mellin convolution formula 
for a ratio, namely 

g(u)= I fa£)(l-l)"-ifa( v )"f( v )dv (1.6) 
Jv<u a a u 

in the scalar variable case, and for the matrix variable case 

g(U)= [ 1 (y5 f /-iyl)_L^/_yl[/-iyl|<*- 
Jv<u l p\ a ) 

x\Vf£\U\-W>f{y)4V. (1.7) 

Definition 1.2. Fractional integral operator of the second kind of order a for the one variable real scalar 
case or one matrix in the real case Let f\ and fi be as defined in Definition 1.1. Then the fractional integral 
operator of the second kind of order a. is defined as the Mellin convolution of a product and given by 

«?(«)= / -M-)^-Al--) a - 1 Mv)f(v)dv (1.8) 

Jv>u v v r («) v 

for the real scalar case and for the real matrix case 

ff(i0 = F7T / \v\- e Z L Mv- i uv-l)\i-v-±uv-i\ a -t L <h(v)f(v)dv. (1.9) 

Note that if fi(xi) and 72(^2) are statistical densities of real positive scalar random variables xi and X2 
then <7(u) will represent the density of the ratio u — In the matrix case g(U) will represent the density of 

U, where X\ = vilf~ 1 Vi and X 2 = V. This is the statistical connection and later we will see that Kober 
operators are constant multiples of statistical densities. 



Special cases 
Case (1): Let 



pi(x) = x[ 1 , (j)2{x 2 ) = 1 then 

^) = frr A-) c - 1 (i--r- 1 4/(^ 

1 (a) J v u u u z 

y—C—a f 

= T7^T / («-") a - 1 « c /(«)d« 

1 \<Ji) J v <u 



Kober operator of the first kind for 5i(a) > (1-10) 



In this case it is easily seen that r ^ a ^ g(u) is the density of u = ^ where x\ and x 2 are statistically 
independently distributed real scalar positive random variables. 

Case (2): Let 

<j>i(x\) = x\ _1 , <p 2 (x 2 ) = x^ c , ( = -a. Then 
^ = TT7 /" (u-v) a - 1 f(v)dv = a D- a f(x) 

1 («) Jv=a 

= J* fix) 

= Ricmann-Liouville left sided fractional integral of order a. (1-H) 



If a = then we have a special case of the Riemann-Liouville left sided integral operator. When a — >• — oc 
then it gives - 00 W~ a f(x) = left sided Weyl operator. 

Case (3): 

4>! (an) = x{~ 1 2 F 1 (a + 0, -7; a; (1 - -)) 

u 

then g(u) gives the Saigo operator of the first kind. 
Case (4): 



0i (zi) = a; J, ^2(3:2) = 1 then 

V - c - a (v-u) a - 1 f{v)dv (1.12) 



= K^' a f(u) = Kober operator of the second kind (1-13) 

This can be interpreted as a statistical density. In fact, ^7^~jT^ ff( M ) * s a statistical density of the product 
u = where x\ and X2 are two positive statistically independent real scalar random variables with 

densities f\(x{) and f 2 (x 2 ) respectively where /1 is a type-1 beta density with the parameters (£ + I, a). 

Case (5): 

(f>(xi) = l,(j) 2 (x 2 ) = x 2 then 

T(a) J v v v T(a) J v>u 

= right sided Weyl fractional integral of order a. (1-14) 

Note that we can have the matrix-variate cases also corresponding to the cases (1) to (5) above. Since they 
are straight forward the details are not given here. 

Definition 1.3. Fractional integral operator of the first kind of orders a\, ...,ak in the multivariable case 
Let fi(xi, Xk) and f 2 (v\, .., Vk) be real- valued scalar functions of the scalar variables x\, xu and v\, Vk 
respectively, where f\ is of the form 



Then the fractional integral operator of the first kind of orders a±, a k in the multivariable scalar case is 
given by 



x f( Vl ,...,v k )dV,dV = d Vl A... Adv k . (1.15) 
In the corresponding many matrix variable case /i is given by 

fi(u u ...,u k ) = Mvfu^v^^vfu^v*) 

j=1 l p{ a 3) Jv 3 <U 3 

x f 2 (y lt V k )dV, dV = dV 1 A ... A dV k . (1.16) 

Definition 1.4. Fractional integral operator of the second kind of orders ai,...,a k in the multivariable 
case Let fi(x\, ...,x k ) and fi{v\, .-,Ufe) be real-valued scalar functions of the scalar variables x\, ...,x k and 
Vi,...,v k respectively, where f\ is of the form 

Then the fractional integral operator of the second kind of orders ai,...,a k in the multivariable scalar case 
is given by 



x —)/(«i, "fc)dV,dV = d«i A ... A dv k . (1.17) 

«1 ^A: 



In the corresponding many matrix variable case <? is given by 



g(m,...,U k ) = f ... f h UiV x i ,...,v k h U k V k *) 

Jyi>c/i Jv k >u k 



x f(V u ..., V k )dV, dV = dVi A ... A dV k . (1.18) 

We will give the following formal definitions for fractional derivatives. Fractional derivatives will be 
defined as the following fractional integrals. To this end let us have the following notations for the fractional 
integrals. All the following integrals are of order a in the single variable case and of orders ot\, ...,a k in the 
multivariable case. Let 



IF" — fractional integral of the first kind, one variable case, eq: (1.6) 

IS" — fractional integral of the second kind, one variable case, equation, cq: (1.; 

I Fx = fractional integral of the first kind, one matrix variable case, eq: (1.7) 

IS X = fractional integral of the second kind, one matrix variable case, eq: (1.9) 
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IFxi,'.'.'.\'xk = fractional integral of the first kind, cq: (1.15) 
ISxl','.'.'.',xu = fractional integral of the second kind, eq: (1.17) 
IFxl' 'x k k = fractional integral of the first kind, eq: (1.16) 

IS x 1 '"'' < x K — fractional integral of the second kind, eq: (1.18) (1-19) 

For all the following definitions m — [a] + 1 where [a] denotes the integer part of 3?(a) > and mj = [ctj] + 1 
where [a^] denotes the integer part of 5i(a!j) > 0, sj = l,...,fc. Fractional integral can be written as an 
antiderivative and then, for example, the fractional integral IF" can be written as an antiderivative as 
IF" — DF~ a . Let -D™ denote the integer order derivative with respect to the real scalar variable x and of 
order m,m — 0, 1, 2, .... Then symbolically we can write 

D™[IF™~ a } = D™[DF-( m - a 1] = DF* (1.20) 

the corresponding fractional derivative of order a. We can take (1.20) to define fractional derivative, where 
let m = [a] + 1, with [a] denoting the integer part of jR(a). Let D™ T = denote the partial derivative of 
order m r with respect to Xj. Then, for example, we can define fractional derivatives of orders ai,...,a k in 
the scalar multivariable case as 



Jjmi jym k jpmi—cti,...,m k — a k _ ]J m i JJ m k j^p-(mi-ai),...,-(m k -a k ) 

1 '*' \z X\,...,X k 1 *'* X\....,X k 

Xl,...,X k 

In the matrix-variate case we will introduce the following operator. Consider the matrix of partial derivatives, 

, d N <) 



i >'-y (i-2i) 



dX= { ^ fel determinant of ^ (1.22) 



and the corresponding matrices and determinants of partial differential operators. Then in the single matrix- 
variate case we can define the fractional derivative of order a as follows: 

| j9_ rji ^-a = \£_ rDF - {m - a) = Dp a (L23) 

If there are several matrices then we consider the dctcrmiants of the matrices of partial differential operators 
as the differential operators and, for example, we can define 



d , m . , d 



I mi 



m k jpm 1 —ai,...,m k -a k 



OX: - i)X L V V - 

_ ® I mi \_^_\m k n n-(mi-ai),...-(m l -a t ) 

"fe 1 - l dx~ k l x v - 

= df^;;;;X (1.24) 

where mj = [ctj] + l, j = l,...,k. If any aj is a positive integer then the corresponding derivative is a 
(nij — l)th order total partial derivative and not a fractional one. 

Definition 1.5. Fractional derivative Let the "/" in (1.19) be replaced by "£)" to denote the corresponding 
fractional derivative. Then the fractional derivative of order a for the single scalar variable case will be defined 
as in (1.20) with the corresponding definitions for all other single variable fractional integrals in (1.19). The 
fractional derivative of order a in the single matrix case will be defined as in (1.23) with the corresponding 
changes for all the ssingle matrix variable cases in (1.19). Fractional derivatives of orders cti, ...,ak in the 
multivariable case will be defined as in (1.21) with the corresponding changes for all cases listed in (1.19), 
and in the many matrix variable case will be defined as in (1.24) with the corresponding changes for all the 
cases in (1.19). 
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2. Fractional Operators in the Many Matrix-variate Cases 

Let X\, Xk and Vi, Vk be two sequences of p x p matrix random variables where between the sets 
the two sets are statistically independently distributed. Further, let Xi,...,Xk be mutually independently 
distributed type-1 real matrix-variate beta random variables with the parameters + ^-,otj),j — 1, k. 
That is, Xj has the density 



r p("i) r p(0 + V) 

for O < Xj < I, > 3?(Cj) > — 1 and fj(Xj) = 0, j = 1, fc elsewhere. Consider the transforma- 

— — — — p+i p+i 

tion V- 2 UjVj 2 , j = 1, ...,k then the Jacobian is Vi | s— . . . | v^, | s— . Substituting in (2.1) the joint density 

of Ui, U k , denoted by g{Ui, U K ), is given by 



m rn r A ^ + 6 + ^) 



i r p (ay)r p (o + 

JVi>Ui 



_ f A r P K- + Q + ^) c , 
7 = L 1 r p (a J )r p (o + £±i) 1 jl 

x/ l^l-^-^l^-^r- 2 * 1 }/^!,...^^. (2.2) 

Hence we will define Kober operator of the second kind and of orders (ai, a^) for the many matrix-variate 
case, and denoted as follows: 



K^^- k f{U u ...,U k ) = {H JMi f Vj « «*|Vj u^-^} 

x f(V u ...,V k )dV. (2.3) 
Therefore this Kober operator is a constant times a statistical density function, namely, 

{ p^tjtL wi. ■•■> ^) - 4?;.^f 1, - , V(^i, (2-4) 

l p (aj + (,j + 2 J 

Now, let us consider f\ having a joint density of the matrix variables X\, Xk and /2 is a joint density of 
the matrix variables V\,...,Vk where the two sets are independently distributed. Then we can have several 
interesting results where the Kober operator of (2.4) will become constant multiples of statistical densities 
coming from various considerations. 

Theorem 2.1. Let the two sets X\, Xk and V\, ...,Vk of matrices be independently distributed. Further, 
let Xi, Xk have a joint type-1 Dirichlet density with the parameters (Q + 2^— , otj), j = 1, k. Consider 
the transformation 

X 1 =Y 1 

X 2 = (I - Y 1 )^Y 2 (I - Yi)i 

X 3 = (I - Y^...(I Yi)^.(J - Yi)*...(J - Kf-i)', j = 2, k. (2.5) 

Or 



6 



Y j = (I-X 1 -...-X j _ 1 )-ix j (I-X 1 -...-X j _ 1 )-i,j = 2,...,k,Y 1 =X 1 . (2.6) 
Consider the transformation 

U i = V } Y i V h Y i = Vf^UjVf^j = !- k - ( 2 - 7 ) 

Then the joint density ofU\, Uk is constant times the generalized Kober operator of the second kind defined 
in (24). 

Proof: Under the transformation in (2.5) or (2.6) the Jacobian is 

J =\I — Yip-iX^)!/ _ y 2 |(*-2)(**i)...|/ _ Yk \^ (2.8) 

and that Y\,...,Y k are independently distributed as type-1 real matrix- variate beta random variables with 
the parameters ((j + ^-,Pj),j = 1, k where 

Pj = O+i + 0+2 + +Ck + (k- j) (2.9) 

see, for example, Mathai (1997). Now, it is equivalent to the situation in (2.3) and (2.4) with Yj's standing 
in place of the independently distributed Xj's and hence from (2.4) we have the following result: 

{ n r^fiYS^ 9 ^ Uk) - x&^-'m, u k ) (2.10) 

j = i L p\Pj ■+" y + 2 - 1 

where flj is given in (2.9). Hence the result. 

We can consider several generalized models belonging to the family of generalized type-1 Dirichlet family 
in the many matrix-variate cases. In all such situations we can derive the generalized Kober operator of 
the second kind in many matrices. We will take one such generalization here and obtain a theorem. Let 
fi(Xi, ...,Xk) of the form 

f 1 (X 1 ,...,X k ) = C\X 1 \^\I-X 1 \ Pl \X 2 \ c '\I-X 1 -X 2 \ pa ... 

x \X k \^\I-X 1 -...-X kl \^\I-X 1 -...-X k \^ P -^ (2.11) 

where O < X\ + ... + Xj < I,j = 1, k, SR(Cj) > — l,j = 1, ...,fc, C is the normalizing constant. Other 
conditions on the parameters will be given later. In this connection we can establish the following theorem. 

Theorem 2.2. Let Xi, ...,X k have a joint density as in (2.11). Consider the transformation as in (2.5) and 
(2.6) with the Uj 's and Vj 's defined as in (2.7). Let the joint density ofUi, U k be denoted as g(U\, Uk)- 
Let 

Sj = 0+i + - + Cfc + ft- + ft+i + ... + Pk,j = 1, k. (2.12) 

Then 

<II r^flt^M ^ Uk) = 4?,!%f U k ) (2.13) 

j=l 1 P\^3 ' 2 + "J/ 

w/iere 5j is defined in (2.12) and the Kober operator in (2-4)- 

Proof: We can show that under the transformation in (2.5) or (2.6) the Yj's arc independently distributed 
as real matrix-variate type-1 beta random variables with the parameters (Q + 5j),j — 1, k. Now the 
result follows from the procedure of the proof in Theorem 2.1. For Kober operators in the one variable case 
see Mathai and Haubold (2008). 
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Note 2.1. Special cases connecting to Riemann-Liouville operator, Weyl operator and Saigo operator, 
corresponding to the ones in section 1.1 for Kober operator of the second kind, can be obtained in a parallel 
manner and hence they will not be repeated here. 

3. Fractional Integral Operators of the First Kind in the Case of Many Matrix Variables 

Let fi(Xi, X k ) be a function of many p xp matrices and / 2 (Vi, V&) be another function of another 
sequence of p x p matrices V\, V k . Let /i be of the form 

4 r ft?n ^ |(, ^ |J -^ |a, ^ } (3 ' 1} 

for O <Xj < J,R(aj) > ^ »((,■) > ^,3 = l,...,fc. Let X s = V?U^Vi,j = l,...,k. Let 



x |J- VfU^V? r'^lVjl^lUjl-^+V} 

x f(V u ...,V k )dV (3.2) 

= I$;$l =1, -" ,k Wu...,U k ). (3.3) 

Then (3.3) will be taken as the definition of fractional integral operator of the first kind of orders a\, ...,a k 
in the many matrix variable case. 

Simplifying (3.2), we have a definition for generalized Kober operator of the first kind in many matrix 
variables case. We note that 



x f(v u ...,v k )dvM^) > ^Md) > P -^- (3-4) 

Hence from (3.1) to (3.4) we can have the following theorem. 

Theorem 3.1. Let Ui,...,U k be independently distributed as real matrix variate type-1 beta matrices with 
parameters (Cj > ) > J = 1, k, $t(ttj) > ,5R(Cj) > LetVi,..,V k be another sequence of matrices 

having a joint density f(V±, V k ). Let the two sets (X\, X k ) and (Vi, V k ) be independently distributed. 
- -l - 

LetXj = Vj 2 Uj Vj 2 , j = 1, k. Let g(U\, U k ) be the joint density of U\, ...,U k . Then the Kober operator 
of the first kind for many matrix variables case as defined in (3.3) is a constant multiple of g(U\, U k ) as 
in (3.2). 

We can also have theorems parallel to the ones in section 2 and the proofs are parallel. Hence we list 
two such theorems here without proofs. 

Theorem 3.2. Let X\, ...,X k have a joint real matrix- variate type-1 Dirichlet density with the parameters 
(Ci, Cfc! Cfc+i)- Consider the the transformation in (2.6) and letYi,...,Y k be defined there. Let V\,...,V k be 
another sequence of matrix random variables having a joint density f(Vi,...,V k ) where let {X\, ...,X k ) and 

(Vi, V k ) be independently distributed. Let Yj — Vj 2 Uj V 2 ,j = 1, k. Let the joint density of U\, U k 
be denoted by g(Ui, U k ). Let 

7j=0+i+- + Cfc+i- (3-5) 
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Then 



k 

(3.6) 



9 {u lt u k) = { n ^^ ^af'"^ - u *) 

k 

HI t¥t^M Ui > Uk) = ^t^f^ (3-7) 



or 



Theorem 3.3. Le£ have a joint density as in (2.11) and the remaining transformations and 

notations remain as in Theorem 3.2. Let 

Sj = Cj+i + ••• + Cfc+i + Pj + ••• + /3 k - (3.8) 
Let the joint density of U\, Uk be again denoted by g(Ui, U k ). Then g{U\, Uk) is a density and 

<h T^C%) }9iUu Uk) = / fe^f 1 '---' fe /(t/i, Uk). (3.8) 



4. M-transforms for the Fractional Integral Operators in the Many Matrix-variate Case 

Here we look at the M-transforms for fractional integral operators of the first and second kind in the 
many matrix-variate case. Consider the first kind operator in (3.4). The M-transform is given by 



f ...J Pl 

JU!>0 Ju k >o 



>m ■ -^...\u k r - 

lUi>o Ju k >o 



l u lt ...,u, 

The X-integral is given by 



x I$?$i~ 1 ''"' k f(U 1 , U k )dU! A ... A dU k 



I P^-^m-^pj - v^-^duj = j \v i \-<>+"-'¥ j Y]>0 \Y 3 r 

Ju 3 >v 3 

x \i + Y j \-^- s ' +p ^dY h Tj = Uj - Vj,Yj = vrhjvr 1 * 



_ p+1 

S ~~T~ 



,+ Si -e+i r p (aj)r p (£±^ + Q - a,) 
r p (£±± + aj + Cj - sj) 



by evaluating the integral by using type-1 real matrix-variate beta integral, for 5?(s) < 5R(C + ^-),^St(aj) > 
Now the V}-integrals give the M-transform of f(Vi, V k ). Hence we can have the following theorem. 

Theorem 4.1. For the fractional integral of the first kind of orders ai,...,ak defined in (3.4) the M- 
transform is given by 

M{IpS&i =1 ''"' k f(Ui,...,U k )} = f*( Sl ,...,Sk) 

Tp(^2 ^ C; — s j) 



x II i 
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2 

f^T^+aj + Q-Sj) 



for K(sj-) < n(Q + 2 t 1 ), K(aj) > ^-,j = 1, k., where f*(si, s k ) is the M=transform of f(Vi, V k ). 
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In a similar manner we can workout the M-transform of fractional integral operator of the second kind 
in the many matrix- variate case. In this context we start with (2.3). The M-tranform is given by 



UjY'-^- / \VA~ Ci ~ ai 
pW) Ju 3 >o Jv 3 >u 3 



p+l 



- P+l 

2 \ a J-—dUj. 



x \Vj - U j \ ai - — dU j }f(V u V k )dV 1 A ... A dV k 

The C/j-integral is given by 

/ \u j \'>- e P+<*\v j -u i \ a '- e P<iu j = \v j \ a '- e P f w^-^i-v-iujvr* 
■ii . \ .ii \ 

_ i _ i 

Put Yj — V 2 UjVj 2 and integrate out by using a matrix- variate type-1 beta integral 

= ,y |a ,-p+A +Sl+Cl Tp(aj)T p (Q + sj) 
Fpiaj+Q+Sj) ■ 

Now the Vj -integrals give the M-transform of /. Hence we have the following theorem. 

Theorem 5.2. For the fractional integral operator of the second kind defined in (2.3) the M-transform is 
given by 

m^M =1 '"' ,k m, u k )} = /•(*,..., s k ) n r p ^ + o + Bj ) 

/or 3i(aj) > E ^ L , 5R(Cj + Sj) > E ^-, j = 1, fc, w/iere /*(si, s fe ) is the M-transform of f(V\,..., V k ). 

Note that for k = 1 the corresponding M-transforms in the one matrix variable case, for p — 1 the 
corresponding Mellin transforms in the k scalar variables case and for p = 1 , k = 1 the corresponding Mellin 
transforms in the one variable case for the Kober operators of the first and second kinds are obtained from 
Theorems 5.1 and 5.2 respectively. 
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